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Abstract 

Strain and rotation fields of dislocations in monolayer graphene have been mapped in a recent 
experiment. These fields are finite everywhere and differ from those given by linear elasticity which 
does not consider rotation explicitly and predicts infinite rotation and strains at the dislocation 
point. A derivative regularization of two-dimensional linear elasticity fixes these shortcomings. The 
theory adds rotation, dislocation and residual strain energies to the usual elastic energy. There are 
two extra material constants that determine the size of the dislocation core and are determined 
from experimental data. These findings are useful for studies of dislocations in graphene and for 
analyzing incipient plasticity in two dimensional crystals. 

PACS numbers: 61.48. Gh,68.65.Pq,64.70.-p 



1 



I. INTRODUCTION 



Advanced imaging techniques with high-resolution transmission electron microscopes 
have produced accurate images of the displacement and strain fields of dislocations up to 0.03 
A.-i^ In graphene and other 2D crystals, images of dislocation cores as defects in the crystal 
lattice and real-time pictures of defect evolution are now being obtained.-"^ Gliding and 
climbing motion of dislocations in graphene has been demonstrated in recent experiments.- 
In the same experiments, the strain and rotation fields (i.e., the symmetric and the an- 
tisymmetric parts of the 2D displacement vector gradient, respectively) near dislocations 
have been obtained.- These findings have surprising elements that are worth examining and 
offer the opportunity to understand plastic deformation in nanoscale materials. At the mi- 
crometer scale, there are effective computational theories of line dislocations that rely on a 
mixture of theory and empirical rules for dislocation interaction and motion.- Precise mea- 
surements and theoretical understanding at the nanometer scale could help bridging the gap 
between scales. On the other hand, controlling electronic properties in graphene through 
strain engineering is a promising concept^ and dislocations and defects could play an impor- 
tant role. For example, special line defects (grain boundaries) may be used to filter electrons 
from different graphene valleys (valleytronic devices).- Experiments have shown that vacan- 
cies and point defects produce paramagnetism^ and appropriate strain fields induce strong 
pseudo-magnetic fields and Landau levels.— 

Dislocations in two-dimensional (2D) linear elasticity are solutions of the Navier equations 
with point forces at the dislocation point. In 3D, they correspond to planar edge dislocations 
that respond to forces supported on the dislocation line.— Both the strain^^ and the rotation 
are inversely proportional to the distance to the dislocation point and are therefore singular 
there." A simple way to eliminate this unphysical behavior is to regularize elasticity near the 
dislocation cores by using the correct lattice structure of the crystal. This produces spatially 
discrete elasticity, an idea that goes back to Frenkel and Kontorova.— Lattice regularization 
leads naturally to the concept of Peierls stress and of dislocation glide as traveling wave 
motion;— see RefJlGl for the case of graphene. Another widely used regularization is the semi- 
continuous Peierls-Nabarro (PN) model and its generalizations that contain the interatomic 
distance as a parameter.^ The original PN model regularizes the dislocation singularity 
at the dislocation line but produces unrealistically large strains at the dislocation core (cf. 
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page 223 in Ref. Il2f ). More recently, some authors have replaced the vertical coordinate y 
(a variable) instead of the interatomic distance (a fixed number) in the PN expressions for 
y = 0.— The resulting formulas,- including a modification that adds an extra free parameter 
to the PN expressions,^^ fit well the experimental strain field in graphene far enough from the 
dislocation point. ^ However the modified PN formulas restore the singularity of the strains 
at the dislocation point, thereby losing the raison d'etre of the PN model. No explanation 
of the rotation caused by dislocations is offered. 

In this paper, we present a consistent modification of 2D continuum elasticity that reg- 
ularizes strains and rotations at the dislocation core and agrees with experimental data in 



graphene.- See Ref. 18 for related work in 3D. The idea is to add three terms to the usual 
strain energy density in continuum elasticity: (i) a term proportional to the square of the 
rotation, (ii) a term proportional to the square of the dislocation density vector, and (iii) 
a background term containing the contribution of the residual stress needed to equilibrate 
the dislocations. Even if rotating the crystal costs little, the term (i) will forbid the infinite 
twists at the dislocation point produced by linear elasticity.— A side effect of this term is 
that the stress and distortion tensors are not symmetric. The dislocation energy (ii) is not 
zero wherever there are dislocations and it produces the moment stress needed to balance 
the antisymmetric part of the stress tensor. The residual stress tensor in (iii) is a sort of bare 
stress and the whole theory can be thought of as a derivative regularization^^ of continuum 
elasticity: the equations of motion contain higher derivative terms that produce a smoother 
dressed stress. The material constants associated to rotations and dislocations yield two 
lengths, and that render finite the rotations and the strains at the dislocation point. 
Over these lengths, rotations and strains differ from the usual expressions of continuum 
elasticity. We find Iq,, Iq, and the associated material constants by comparing the formulas 
provided by the theory with experimental data.- These modifications of linear elasticity in 
presence of dislocations may pave the way to bridging the nano and micro scales. 
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II. STRAIN ENERGY AND EQUATIONS OF MOTION 



For 2D planar graphene with displacement vector Ui {ui = u, U2 = v), the potential 
energy is split in elastic, rotation, dislocation and background parts 



-2crifcAfe] dx dy, (1) 



where A and /i are the 2D Lame moduli, 7 is the rotation modulus, /3(ifc) = {f3ik + Pki)/'^ and 
/3[i2] = (/3i2 — /32i)/2 are the elastic strain and the rotation, respectively. Sum over repeated 
indices is implied. The displacement gradients, Ui^ = djUi = (3ij + , contain elastic and 
plastic distortions, (3ij and Pfj, respectively. The dislocation strain energy (with associated 
modulus > 0) is quadratic in the dislocation density vector a^. The latter can be extracted 
from the definition of the Burgers vector,— 

— bi = ^ dui = <j) dxjUij = tkj j dkf^ijdx dy, (2) 

where we have replaced Pij = Uij, the line integral is calculated on a contour encircling 
the dislocation point, we have used the Stokes theorem, and e^j is the 2D completely an- 
tisymmetric tensor with eu = 1. From this expression, the dislocation density vector- 
is 

«i = £jkf3ij,k- (3) 

For a point dislocation with Burgers vector bi located aX x = y = 0, the dislocation density 
is just ai = bi6{x)6{y). While ejkUijk = 0, neither the elastic or the plastic distortions are 
gradients of functions. In fact, ([3]) indicates that the incompatibility of the elastic distortion 
equals the dislocation density. The last term in ([T]) (the residual or background energy) 
contains the residual stress tensor determined by the equilibrium condition 

= 0- (4) 

The potential energy is a functional of the distortion tensor In equilibrium, its 
variation satisfies 

6W dW ^ dW ^ 

- — -, r = 0, (5) 



6/3ij (x, y) d/3ij (x, y) d/3ij^k {x, y) 
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which yields 

CTij + <^ekjai,k = CTij, (6) 
where the stress arises from the first three terms in ([T]) 

aij = XPnnSij + "^fiP^ij) + 27/3[jj] 

= 2/i f ^5,, + +27%]. (7) 



Here z/ = A/(A + 2/i) is the 2D Poisson ratio. This imphes 

~ 2^^{l + v) ^^''^ ~ 2/i V""^*'^ 1 + 1/ 

The residual stresses satisfy 0"° ^ = 0. Then differentiating ([6]) and using the compatibility 
identity, ekjai^kj = 0, we get the usual condition of force equilibrium, 

a,,-, = 0. (9) 



III. STRESS, STRAIN AND ROTATION 

To solve (E]), we write the dislocation density functional of the stress and 

calculate the residual stress for a given dislocation configuration. 

A. Equations of motion for the stress 

Equation ([3]) gives the dislocation density in terms of the distortion, which can be calcu- 
lated as a function of stress using ([8]). After some algebra and simplifications through the 



use of aij J = 0, we obtain 



1/1 Ma 

2 V7 /^/ 

Aa 

€kiOli,k = "2,1 - "1,2 = -7^T'^ V' (10) 

2/i(l — U) 



where A/ = VV = {df + dl)f. Using 



^k[j(^i].k = \ ^ ^k{jOti)^k 

0, i= i 



"1,2, 2 = j = 1, 

"2,1, ^=J=2, (11) 

l(ai,i - 02,2), ^T^j, 



that follows from (|3]), the antisymmetric part of ([6]) becomes 

^[12] = (^[12]- 

Similarly, taking the trace in (|6]) and using ( ITO|) and ( ITT]) , we obtain 



4 V/i 7 



(12) 



2/i(l + 



-A 



(13) 



B. Bare and dressed stress functions 

As indicated in Appendix |Al Equations ( !T2|) and ( !T3|) can be solved using the Airy stress 
functions, $ and such that 



CTll = $,22 - ^,12, 0'12 = -^",12 + ^,1 



0"21 



,11' 



-$,12 - ^,22, 0^22 = $,11 + ^,12- 



(14) 



A similar expression holds for the residual stress components, a^p in terms of residual stress 
functions, $° and Then we obtain ann = A$, a[i2] = |A\1/, which, inserted into ( !T2|) 
and (|T3|) . yield the equations 



(l-/|A)$ = <l>o, (l-4A)^ = ^o, 
/2 - ^ /2 _ i ?1 , 1 



(15) 
(16) 



Eq. ( IT5l) indicates that $ and \1/ are higher derivative regularizations of the usual Airy stress 
functions. To solve ( [T5|) for a single dislocation located at the origin, we need to first find the 
residual stress functions $° and According to Eq. ([3]), the residual distortion induced 
by the dislocation satisfies 



^jk(3ij,k = -biS{x)6{y), i.e., 
/3a,2-A°,i = -M(x)%). 



(17) 



Using (|8]) to rewrite (fT7|) in terms of the residual stresses and stress functions, we obtain 



A2$° = -2/i(l + u){hd2 - hd,)5{x)5{y), 



/i + 7 



(Ml+M2)5(x)5(|/). 



(18) 
(19) 



6 



Thus the residual stress functions are 



$0 
^0 



-2/i(l + z/)(6i92-629i)G(x,i/), 
47/i 



/i + 7 



{hidi + h2d2)G{x,y), 



(20) 



G{x,y) = —rHnr, A^G{x,y) = 6{x)6{y). 
an 



$° is the stress function of an edge dislocation in 2D linear elasticity.— The corresponding 
bare stress, strain and rotation are singular at the origin. For a dislocation with Burgers 
vector (&, 0), the solutions of the inhomogeneous Helmholtz equations f|T5l) are 



A 



2 ^ 



dA In r + Al 



In r + 



Inr + Koi — 

tvjr 



(21) 
(22) 



where A = fib{l + u)/ (27r), B = ji'^h/[n{jj, + 7)]. These dressed stress functions become the 
residual (bare) stress functions ( !20l) for lq> = Iq, = Q and coincide with those of a 3D straight 
edge dislocation^ except for the different value of the prefactor A in ( 1211) . From the stress 
functions, we can find the uniformly bounded stress, strain, rotation and dislocation density 
vector. We get 



/3[12] 



h{l + u)y 



27rr2 
fibx 



7r(yU + 7)r2 
2(?/2 — vr'^)r 



1 K \ 

1 ^ K f ^ 



4/^ 



1 - 2i^y + 2x' + ^(r - 3a; 



2{y^ - 3x^)K2 



X -y 



4/2 
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~2^ 



/3(12) 



All 

,2 /^2 q„,2n 



2y^r 



2/ir'* [ " r 
Kj^]+2{x^-3y^)K2 



+B 



1^22 — Pnn — Al- 



(23) 
(24) 



(25) 



(26) 
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2^ 



4 ^ 



Inr + /|AKo 



Inr + llAKo 



-bil + lldl + lldl)5{x)5iy) 



-b{l + lid'y + lidMx)S{y)- 



b 
b 



+ lldlAKo 



iTcr r 



2' 



3(r2-4/|)(x2-y2) 



4r3/| 
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4/| + r 
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3/, 
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3(r^-4/|: 



6L 



(27) 



(28) 



We can calculate the Burgers vector distribution by integrating the dislocation density on a 
disk of radius r centered at the origin. The resulting vector is (&(r), 0) with 



b{r) 




2tt 



b + 



^0 
b 



air drd9 = b+- rdr ( 9^ + -dr 

r 



U 



+ K2 



2 

+ ^0 



+ K2 



+ 



2/<j) 
r 



3r 
2/<j, 



r 



2/,t 3r 



2/, 



+ — -TT^M^l^ -;T^^3^ -TT^i^S^ 



2L 



U 



2L 



(29) 



The distribution b{r) monotonically increases from to 6 as r —t- 00. 



IV. RESULTS 

Figure [1] shows the contours of the strain f3u{x,y) and the rotation f3[i2]{x,y). Fig. [2] 
compares the rotation along the horizontal axis given by f l24|) and the rotation given by linear 
elasticity, which is unbounded as x — )■ 0. We have found the moduli 7 and <^ (equivalently, 
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FIG. 1: Contours of (a) strain (b) rotation /3[i2]- 
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FIG. 2: (Color online) Rotation profiles at y = 0. Solid line corresponds to ()24p and dashed line 
to continuum elasticity, /3[i2] = 6cos 0/(27rr). 



/$ and by fit to experiments. § The strain e^x inspired by the PN model and Ref. [iTI is^ 



^xx = —h{l — u) sinO 

4(1 - uf cos^ e + (2a/ - l)aj sin^ 9 



(30) 



7rr[4(l-z/)2cos2^ + aJsin2^]2 ' 

when written in polar coordinates. See also Supplementary Eq. (2) in Ref. U Note that 
the strain ( 130|) is singular at r = 0, thereby defeating the purpose of the PN model. The 
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best fit to graphene data at different nonzero distances r is reported to be aj = 1.5.- The 
contours of this strain and those of Pu given by fl25|) should agree. On the other hand, the 
rotation about a dislocation (heptagon-pentagon defect) in graphene goes from -0.5 to 0.5 
radians.^ Using 0.5 radians for the maximum rotation in fl2^ fixes the scale We find that 
<j/ (A -|- 2/i) = 6 X lO^^fe^ and 7 = 12/i. The Burgers vector of a dislocation in graphene is the 
lattice constant, b = 1.42^3 A. At 300 K, we have /i = 9.95 eV/A^, A + 2/i = 22.47 eV/A^-^i 
and therefore <^ = 0.08155 eV and 7 = 119.4 eV/A^. Furthermore the characteristic lengths 

= 0.0256 and = 0.0196 are below 0.1 A. The dressed strain and rotation behave similarly 
to the bare ones a short distance away from the dislocation point. At larger distances, the 
empirical formula f l30p with a/ = 1.5 used to fit our theory gives similar strain contours. 

We can compare the results given by the present theory with those produced by periodized 
discrete elasticity. For this purpose, we have approximated the partial derivatives in the 
definition of the distortion tensor by finite differences using the honeycomb lattice. Very few 
values are available to reconstruct partial derivatives from finite differences at the core of 
the dislocation which would produce a pixelated figure (not shown). The resulting density 
plot is similar to that yielded by our regularized continuum model, with positive (negative) 
values of rotation for a; > (x < 0). In the case of periodized discrete elasticity, the 
symmetry y — —y of Fig. [Hb) of the continuum theory is broken due to discreteness effects 
which are very strong near the dislocation point. This symmetry breakdown is also visible 
in experiments; see Supplementary Figures SIO, Sll and S15 in Reffl. 



V. CONCLUSIONS 



We have calculated the strain and rotation contours about a single dislocation in graphene 
by regularizing continuum elasticity in two ways. A lattice regularization (periodized dis- 
crete elasticity)^ eliminates the singularities at the dislocation point and gives contour 
plots that agree with experimental data.- We can also carry out a derivative regularization 
of elasticity by means of a continuum theory that introduces rotation, dislocation and back- 
ground (residual) strain densities. This theory contains two additional material constants 
that we have estimated by using data from experiments. While we have considered a pla- 
nar graphene sheet, off-plane vertical displacements can be incorporated to our continuum 



theory by including a bending energy and using ideas similar to those in Ref. 
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Appendix A: Proof of the existence of the Airy stress functions 



The arguments in this Appendix are adapted from Ref. l23|. Since Uij = f3ij + ([3]) 
imphes that = — where a° is defined by replacing the residual distortion instead 
of the distortion in (IHl) . In (161) , let us denote 



(^tj = ^^kjai,k, (yf^ = ^ekjalk =^ (a^ + crfj) - ((t° + af^) = 0. (A.l) 



Clearly, 



and we can introduce eight functions 0, 0, ip, ip, 4'^, 0°, "0°; "0°; 



an = 4>y, CTu = -(f)x, 0-21 = -02/, 0-22 = 0x, (A. 3) 

afi = tjjy, crfa = -^x, = -i'y, (7^2 = i'x, (A.4) 

with similar definition for the residual stresses. In this Appendix, subscripts in the stress 
functions imply partial derivatives: (p^ = dx4>, etc. To show that these four functions 
can be written in terms of only four stress functions, we proceed as follows. Firstly, we 
show that ip cind ip are proportional to the components of the dislocation density vector. 
i/jy = crfi = — ^ai,2 and —ipx = crfi = imply that ip = —(^ai. Similarly, we can show 

that -0 = ^0:2, "0° = <J«2) = — <J«i- 

Now the antisymmetric part of can be written as cr[i2] + |^efc[2tti],A: — (cr['i2] + |'?ffc[2tti] k) ~ 
0, i.e. 

(^12 - (^21 + ^^("1,1 + "2,2) - o-?2 + 0-^1 - + "2,2) = 0- (A.5) 
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Using (lA.Sp . this equation becomes 

= dy(4>+ 1^2) - 5, (0 - l^ai) - dy (0° + + - l^a?) . (A.6) 

From which we can set 

0+1^^2 = $.., (P-^<;a^ = %, 0O + l^aO = $o, / - lc;a? = ^ 

= $.-|«2, = + 0o = $o-ia°, / = $o + ia?, (A.7) 

in terms of new functions $ and The strain compatibihty conditions UjC^ij = 0, ey-a^j = 
0, imply that the dislocation density vector derives from a potential, and we have qai = —"^x, 
= ~'^y similar expressions for Then (1A.7P becomes 

= $^. + ^^, = $^-iI^^, V = *x, ^ = -^?;, (A.8) 

with similar expressions for the residual stress functions. Substituting this into (1A.3P and 
flA.4p . we obtain 



0-11 = ^yy - "^xy, 0-12 = -^xy + "^xx, ^21 = -^xy " "^yy, (^22 = ^xx + ^xy,(A.9) 

crfi = ^xy, o-f2 = -^xx, 0-21 = ^yy, ^22 = -"^xy, (A.IO) 

with similar expressions for the residual stress components and residual stress functions. 
Equation (1A.9P is the same as ffUl) . If = 0, (1A.9P is the usual relation between the stresses 
and the Airy stress function of linear elasticity. 
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